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Measurements of the magnetic irreversibility line and 
time-logarithmic decay of the magnetization are described for 
three FeiOz samples composed of regular amorphous, acicu- 
lar amorphous and crystalline nanoparticles. The relaxation 
rate is the largest and the irreversibility temperature is the 
lowest for the regular amorphous nanoparticles. The crys- 
talline material exhibits the lowest relaxation rate and the 
largest irreversibility temperature. We develop a phenomeno- 
logical model to explain the details of the experimental re- 
sults. The main new aspect of the model is the dependence 
of the barrier for magnetic relaxation on the instantaneous 
magnetization and therefore on time. The time dependent 
barrier yields a natural explanation to the time-logarithmic 
decay of the magnetization. Interactions between particles as 
well as shape and crystalline magnetic anisotropies define a 
new energy scale that controls the magnetic irreversibility. In- 
troducing this energy scale yields a self-consistent explanation 
of the experimental data. 



I. INTRODUCTION 

A ferromagnetic particle becomes monodomain if its 
size d is reduced below a critical value d cr ~ 1 — 100 nm, 
determined by the competition between dipole and ex- 
change energiesEm. Below this critical size, the energy 
loss due to creation of magnetic domain walls (propor- 
tional to d 2 ) is larger than the gain due to disappearance 
of the dipole magnetic field energy (proportional to d 3 ). 
Such monodomain ferromagnetic particles can be viewed 
as large magnetic units, each having a magnetic moment 
of thousands Bohr magnetons. Usually neighboring par- 
ticles are well separated (10 — 20 nm), and direct ex- 
change between particles may be neglected. Thus, the 
magnetic properties of an assembly of nanoparticles are 
determined by the dipole field energy alone, w^th thermal 
and magnetic anisotropy energies (see e.gjjia). 

Experiments conducted on magnetic nanoparticles 
show irreversible magnetic behavior below the "irre- 
versibility line" Ti rr {H). In particular, the zero- field 



cooled (ZFC) and field cooled (FC) magnetization curves 
do not coincide, and magnetic hysteresis appears in M 
vs. H curves (see e. gJ3~E3). Moreover, time-logarithmic 
magnetic relaxation, towards the thermodynamic equi- 
librium state, is observed below Ti rr (H). Similar ob- 
servations are reported here for three systems of Fe^O-i 
nanoparticles with different shape and crystalline mag- 
netic anisotropies. These nanoparticle samples were pre- 
pared by a sonochemical meth^cL-which produces " regu- 
lar" amorphous nanoparticlcstH Sonochemical irradi- 
ation carried out in the presence of magnetic field-cesults 
in synthesis of acicular amorphous nanoparticlestJ. An- 
nealing of amorphous particles leads to crystallization. In 
this manner we have prepared regular amorphous, acic- 
ular amorphous and crystalline ferromagnetic nanopar- 
ticles. This enables a study of the effect of shape and 
crystalline anisotropies on magnetic irreversibility and 
the relaxation rate. We find that, qualitatively, all three 
samples exhibit similar irreversible magnetic behavior. 
However, their irreversibility lines and relaxation rates 
differ significantly. 

Irreversible magnetic behavior similar to that de- 
scribed here is also observed in other systems. 
A noticeable example is magnetic irreversibility in 
superconductorsE-3 E3. In such systems the origin for ir- 
reversibility is the interplay between thermal energy and 
some energy barrier, which prevents magnetic reorgani- 
zation in those materials. The microscopic origin of the 
barrier, however, depends on the system. The magnetic 
irreversibility in nanoparticles is conventionally associ- 
ated with the energy required for a particle moment re- 
orientation, overcoming a barrier due to magnetic shape 
or crystalline anisotropy. It is important to note that the 
barrier is considepejd|-tp be independent of the magnetic 
moment itselfail3'tIrE2l. In superconductors, magnetic ir- 
reversibility is due to the inevitable spatial fluctuation of 
the superconducting order parameter caused by defects, 
imperfections etc.; the barrier is the energy required to 
overcome the pinning due to this disorder. 

An important concept in the theory of irreversible 
magnetic properties xtf superconductors, based on the 
work of AndersorOllj, is that the effective barrier for 
magnetic relaxation increases with time. This is because 
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the superconducting shielding current (proportional to 
the magnetization) decays with time, causing a decrease 
in the Lorentz force which drives the fluxons sway of 
their positions. In a nice paper Lottis et alE3 have 
put forward similar arguments to study slow dynamics. 
They noticed the close analogy between ferromagnetic 
assemblies and superconductors. Analyzing the results 
of numerical computations they concluded that the de- 
cay of the demagnetizing field is the origin of what the 
"quasi-logarithmic" relaxation. Although the distribu- 
tion of particle sizes may explain quasi-time-logarithmic 
relaxation in a limited time-interval, it is not necessary 
for the explanation of the experimentally observed time- 
logarithmic relaxation. This approach was later em- 
ployed in other studies of the magnetic relaxation, for 
example, in thin magnetic filmso. 

In this work we adopt the concept of time varying 
barrier and derive the phenomenological model to ex- 
plain magnetic irreversibility and logarithmic magnetic 
relaxation in nanoparticlcs. The physics for the time 
dependence of the barrier in nanoparticles is related to 
the fact that the effective barrier for reorientation of 
the magnetic moment of each nanoparticle depends on 
the internal magnetic field, which includes the average 
dipole field from surrounding nanoparticles. This aver- 
age dipole field decreases with time due to the increase of 
randomness in the orientation of the magnetic moments 
of the surrounding nanoparticles. This, in turn, causes 
the increase of the effective barrier with time, yielding 
a natural explanation to the experimental observation 
of time-logarithmic relaxation, and a sample dependent 
irreversibility line and relaxation rate. Interactions be- 
tween particles as well as shape and crystalline magnetic 
anisotropics define a new energy scale that controls the 
magnetic irreversibility. Introducing this energy scale 
yields a self-consistent explanation of the experimental 
data. 

This article is organized in the following way. In Sec- 
tion II we describe the preparation of the three nanopar- 
ticle systems. We then describe our experimental results 
of irreversible magnetic properties at various tempera- 
tures, fields and times. In Section III we describe our 
phenomenological model and derive equations for the ir- 
reversibility line and the magnetic relaxation. In Section 
IV we compare the predictions of our model with the 
experimental results. 

II. EXPERIMENTAL 

A. Sample preparation and characterization 

Three Fe^O^, samples composed of regular amorphous, 
acicular amorphous and crystallinCj-jripjoparticles were 
prepared by a sonochemical methodlinl3. For the ultra- 
sound irradiation we used VC — 600 Sonics and Materials 
sonicator with Ti horn at 20 kHz and 100 W ■ cm~ 2 . In 



Table I we summarize their features. One molar solution 
of Fe{CO)$ in decaline was sonochemically irradiated for 
three hours in ambient pressure at °C. The powder 
obtained was centrifuged, washed repeatedly with dry 
pentane (6 — 7 times, 8500 rpm), and dried in vacuum 
at room temperature for three hours. The material ob- 
tained has been accumulated from 2 — 3 sonications and 
the total amount of Fe20% was mixed to ensure the re- 
liability of the results. Then, in order to remove organic 
residue, material was annealed in vacuum at 140 — 150 °C 
for three hours. Heating up to this temperature was nec- 
essary to evaporate residua of solvents, particularly de- 
caline which has high a boiling point (189 — 191 °C). The 
amorphous nature of the particles is confirmed by X-ray 
diffraction, differential scanning calorimetry (DSC) anal- 
ysis and electron-diffraction patterns at selected areas as 
shown in Figs, la, 2a, and the inset to 3a, respectively. 
The absence of Bragg peaks in Fig. la demonstrates 
the absence of the long-range order in the atomic struc- 
ture; the large endothermic peak in Fig. 2a indicates an 
amorphous to crystalline transition at ~ 400 °C. The 
electron diffraction pattern of the inset to Fig. 3a also 
confirm the amorphous nature of the particles. A typical 
particle size of ~ 50 nm is inferred from the transmission 
electron micrography (TEM) picture of Fig. 3a. 

Acicular amorphous particles have been prepared by 
performing sonication in external magnetic field of 7 kG 
for three hours. The sonication has been carried out 
in the 0.25M solution of Fe(CO)s in a flask open to 
air. We then repeat the wash and dry procedure as de- 
scribed above. The amorphous nature of the particles 
was confirmed by X-ray diffraction, DSC, and selected 
area electron-diffraction patterns as shown in Figs. 16, 
26, and the inset to 36, respectively. A typical particle 
length of ~ 50 nm and diameter of ~ 5 nm are inferred 
from TEM picture of Fig. 36. 

Heating of amorphous Fe%0§ up to 370 — 380 °C in 
ambient atmosphere for 3 — 4 hours resulted in crystalline 
7 — Fe20s nanoparticles. The nature and the internal 
structure of the crystalline iron oxide were determined 
using the X-ray diffraction shown in Fig. lc. The DSC 
data, Fig. 2c, do not show any endothermic peak. The 
TEM image of Fig. 3c show particles of mean size of 
~ 200 nm. 

The second column in Table I summarizes typical par- 
ticle size for the three samples. The third column in- 
cludes the total surface area of the particles, as measured 
by BET absorption using N2 gas as absorbent. 

B. Magnetic measurements procedure 

A Quantum Design MPMS SQUID magnetometer was 
used for all magnetic measurements reported here. The 
irreversibility line was determined from ZFC and FC 
magnetization measurements. Before taking a data point 
temperature was stabilized with 0.05 K accuracy and a 
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30 sec pause was sustained. The temperature at which 
ZFC and FC merge for a constant field H is defined as 
irreversibility temperature Tj rT . (H) . We define the merg- 
ing point using a criterion \Mpc — Mzfc\ ~ 0.1 emu/g. 

The procedure for measurements of magnetic relax- 
ation at different temperatures is as follows: The sample 
is cooled in H = 2 Tesla from a room temperature (larger 
than Ti rr (2 Tesla)) to a target temperature T, the mag- 
netic field is then reduced to 500 G and the magnetic 
moment is measured for approximately two hours. The 
first data point is taken approximately two minutes after 
the field change. 

The field dependence of the magnetic relaxation rate is 
measured at T = 20 K. At this temperature the field is 
ramped up to H = 2 Tesla and reduced back to a target 
field H, from where the measurements start. The same 
has been repeated for negative field H = — 2 Tesla with 
consequent a increase of the magnetic field to a target 
value and measurements of the magnetic relaxation. 

C. Results 

The experimental results in this section are organized 
as follows: we first show M (T,H = constant) data, and 
related measurements of magnetic relaxation at different 
temperatures. From the merging point of the ZFC and 
FC magnetization curves we extract the irreversibility 
line for the three samples. From the relaxation measure- 
ments we deduce the relaxation rate, as a function of tem- 
perature, for the three samples. We then present mea- 
surements of magnetization loops M (T = constant, H) 
and magnetic relaxation at different values of external 
field. The relaxation rate, as a function of field, is then 
deduced for the three samples. 

Fig. 4 exhibits typical results of ZFC-FC magnetiza- 
tion curves and magnetic relaxation at 500 G for the sam- 
ple composed of amorphous round nanoparticlcs. The 
vertical lines of open circles in Fig. 4 depict the relax- 
ation measurements at different temperatures. The ver- 
tical arrow indicates the direction of the time increase. 
The magnetic moment relaxes towards the equilibrium 
moment M revi determined by the FC curve. In the inset 
to Fig. 4 we zoom out at the ZFC-FC curves and indicate 
by an arrow, the experimental definition of Ti rr . 

The magnetic relaxation data of Fig. 4 are re-plotted 
in Fig. 5 as a function of time. The solid lines in Fig. 
5 are linear fits for M oc In (t) . A qualitatively simi- 
lar time-logarithmic decay is also observed in the other 
two samples. Quantitative differences will be discussed 
below. 

We define the "normalized relaxation rate" R = 
\dM/d\n{t) \ /M c , i. e., the logarithmic slope of the re- 
laxation curve normalized by the magnitude of the ir- 
reversible magnetization at which the relaxation starts, 
M c = Mq — M rev . Here Mq is the initial value of the 
total magnetic moment and M rev is the magnetic mo- 



ment corresponding to a field cooling in 500 G. Fig. 6 
summarizes the values of R as a function of temperature, 
for the three samples. At low enough temperatures, R 
is the lowest for the crystalline sample, intermediate for 
the acicular amorphous sample and the largest for the 
regular amorphous sample. Note, that at higher temper- 
atures it looks as if R (T) curves will cross. This is due 
to a large difference in the absolute values of T irr (90 K, 
162 K and 216 K at 500 G for regular amorphous, acic- 
ular amorphous and crystalline, respectively). As shown 
in Fig. 16, M (T) curves scale with T irr and, therefore, 
in the inset to Fig. 6 we plot R vs. T/Ti rr . In this pre- 
sentation, the whole R(T/Ti rr ) curve of the crystalline 
sample is lower than that of the acicular amorphous sam- 
ple and both are lower than the R curve of the regular 
amorphous sample. 

In Fig. 7 we compare the irreversibility lines for the 
three samples. The largest irreversibility is found in a 
crystalline sample, intermediate in the sample with aci- 
cular particles and the lowest in the regular amorphous 
sample. We explain these observations below. 

Magnetic irreversibility below T irr is also demon- 
strated by measuring the magnetization loops M (H). As 
an example, we show in Fig. 8 M (H) for the amorphous 
nanoparticles at T = 5 and 100 K. Magnetic hystere- 
sis is apparent at 5 K, whereas the behavior is purely 
reversible at 100 K. 

The relaxation at different values of the external mag- 
netic field is shown in Fig. 9. The vertical lines represent 
M (t) curves shown along with the standard magnetiza- 
tion loop. The field dependencies of the relaxation rates 
for our samples are shown in Fig. 10. There is an appar- 
ent change in R between low and high fields. At lower 
fields R is the largest in an amorphous sample, whereas 
at large fields the relaxation rate in an amorphous sample 
is the lowest. 



III. MAGNETIC RELAXATION IN THE 
ASSEMBLY OF NANOPARTICLES 

A. Time dependent effective barrier for magnetic 
reorganization 

Magnetic relaxation is a distinct feature of systems 
with interacting particles, far from thermodynamic equi- 
librium. In an assembly of ferromagnetic nanoparticlcs, 
the elementary process of a change in the magnetization 
is the rotation of the magnetic moment of a nanoparticle 
(or cluster of such magnetic moments). In the following 
we assume that the magnetic anisotropy of each nanopar- 
ticle is strong enough to utilize an Ising-like model, i. 
e., the magnetic moment of each particle is aligned only 
along the anisotropy axis. In Fig. 11a we illustrate 
schematically the orientation of the elementary magnetic 
moments of several of such nanoparticlcs. The full ar- 
rows represent the size and direction of each magnetic 
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moment. The experimentally measured magnetic mo- 
ment is determined by the sum of the projections of each 
individual particle's moment on the direction of the ex- 
ternal magnetic field. Note that the directions of the 
easy axes are randomly distributed. For such a system, 
the energy W of each magnetic nanoparticle, neglecting 
for the moment the interparticle interactions, varies with 
the angle asEx 

W = KV sin 2 [ip-ff]- M p H cos (<p) (1) 

Here is the angle between the easy axis i? and the ex- 
ternal magnetic field I?, and tp is the angle between the 
particle magnetic moment M p and the external field. In 
order to have any magnetic irreversibility and relaxation, 
the KV term in Eq. |] must be larger than the M P H term 
and we will consider this limiting case. The reduced en- 
ergy W/KV of Eq. [l| is plotted in Fig. 12 as a function 
of the angle ip for two different fields Hi = 2.5KV/M P 
(bold) and H 2 — 0.5KV/M„ (light). Since the magnetic 
anisotropy has no preferable direction, there are two min- 
ima in the angular dependence of the energy, as shown in 
Fig. 12. The external magnetic field fixes the direction 
of the lowest minima. We denote by U12 the barrier for 
reorientation from the lowest minima (Wi) to the other 
minima (W2). The backward reorientation requires over- 
coming the energy barrier JJ% \ . 

In order to take inter-particle interactions into account 
we view the field H in Eq. |l] as the internal magnetic 
field, which is the sum of the external field and the dipolc 
field from the surrounding nanoparticles. This local mag- 
netic field depends, on the directions of neighboring mag- 
netic momcntsOcJ. Since the magnetic moment is a sta- 
tistical average of those moments, the local field depends, 
on the average, on the total magnetic moment. This in- 
duces a feedback mechanism: each reorientation of an 
individual nanoparticle decreases the total magnetic mo- 
ment. This is illustrated in Figs. 11a. and 116. Fig- 
ure 11a represents a snapshot of a field-cooled system 
of nanoparticles in which most of the individual mag- 
netic moments are favorably oriented in a direction such 
that their projections are along the external field. Af- 
ter a field decrease, as a result of thermal fluctuations, 
some magnetic moments reorient so that their projection 
is anti-parallel to the external field. The open arrows in 
Fig. 116 represent those reoriented moments. 

Since the local dipole field decreases during this pro- 
cess, the average barrier U12 increases. As indicated in 
the Introduction, an increase of the barrier with time 
is a characteristic of other irreversible systems, such as 
type-II superconductors in the process of magnetic flux 
creep. 

The dynamics resulting from such a scenario is 
sketched Fig. 13. Immediately after reducing the mag- 
netic field, individual magnetic moments arc still along 
the direction of the external field, i. e., in minima Wi of 
Fig. 13, as depicted by the population of the black dots. 
During the relaxation process magnetic moments flip to 



the minimum W2 in the figure. Since, as discussed above, 
this barrier depends on the total magnetic moment via 
dipole fields, it will increase with time as shown in the 
figure, with dipole fields working on the average against 
the external field. The total magnetic moment along the 
magnetic field is thus decreased, as sketched in Fig. 116. 

B. Equations of magnetic relaxation 

In a realistic sample, the directions of easy axes are 
randomly distributed, the particles cannot physically ro- 
tate (e. g., in a dense powder of ferromagnetic nanopar- 
ticles) , and dipole interactions are strong. We will model 
this situation as outlined below. 

Any given particle i in Fig. 11 has an anisotropy axis 
at a fixed angle Oi relative to the external magnetic field. 
The magnetic moment of this particle is then oriented 
at an angle ipi to the field. This angle is defined by 
the non-local energy minimization, due to dipole fields 
of the surrounding. It is important to note that each 
particle interacts with a local magnetic field Hi which 
is the result of a vector sum of the external and dipole 
fields. At small enough external field and large enough 
anisotropy tp+ may have two values: ipi ~ Oi or pi « 
Oi + 7T, which leads to the situation described in Fig. 12, 
with two energy minima at WI « —M p HiCOs(0i) and 
W$ ~ M p Hi cos (Oi). Thermal fluctuations may force 
particle moment in the minima WJ to change its direction 
to another minima W 2 and vise versa. The W\ — » W 2 
rotation requires overcoming a barrier U\ 2l and a barrier 
U{ 2 f° r backward rotation, see Eqs. [l6] and |lj of the Ap- 
pendix, respectively. We then assume that the field Hi 
can be represented as a simple sum of the external field 
H and the collinear to H dipole field Hd (i.e. indepen- 
dent of Oi). The amplitude of a dipole field Hd at any 
given site depends upon orientations of the moments of 
the surrounding particles. If those orientations are to- 
tally random (minima Wi and W 2 are equally occupied) 
the dipole field is small, whereas if all surrounding par- 
ticles are in one of the minima the resulting dipole field 
is maximal. From this simple analysis, we conclude that 
the magnitude of a dipole field depends upon the total 
magnetic moment M. 

Considering the balance of forward and backward ro- 
tations, and averaging over the volume of the sample, 
we show in the Appendix that magnetic relaxation is de- 
scribed by a differential, equation similar to that derived 
for superconductoral3il3: 

where A is an attempt frequency and U is an effective 
barrier for magnetic relaxation given by 

u = u °( i -w) (3) 
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where 



U = 2KV + AM P (H - ~/M rev ) /it 



(4) 



and M = i ( 



ttKV 

2M„ 



Here 7 is the con- 



stant accounting for the strength of the dipole-dipole in- 
teractions, M p is the magnetic moment of an individual 
particle, K is the anisotropy constant and V is the par- 
ticle volume. Apparently, as 7 — ► the energy barrier 
U — * Uq, thus Uo is the barrier in the assembly of non- 
interacting particles. It is worth noting that in our model 
the barrier U depends on the magneticj_cnoment in the 
same way as that used by AndersonliTll3 for a descrip- 
tion of magnetic relaxation in superconductors. 

In the following we analyze magnetic relaxation de- 
scribed by those equations. 

If the barrier for a particle moment reorientation does 
not depend on the total magnetic moment, i.e., 7 = 
and U — Uq, direct integration of Eq. @ yields: 



M = M c exp(-</r) 



(5) 



where M c is the initial irreversible magnetization and 
t — exp(£/o/T) I A is the macroscopic characteristic re- 
laxation time. This result is very similar to that derived 
in early works forj-classical Neel's superparamagnetic re- 
laxation, see e.g., EH. This exponential decay is observed 
experimentally, for example in the work of Wegrowe et 
al.cB on a single nano-wire. 

If interactions are not negligible, Eq. |^ may be rewrit- 
ten in dimcnsionlcss form: 



du 



■ exp (— tt) 



where u — U /T and f = t/t with 
M T 1 ttT 



t 



T 



M c U A 4rfAM p M c A6 



(6) 



(7) 



where we introduced a new energy scale 0, which, as we 
show below, determines the relaxation process and the 
irreversibility line: 



9 = MJJ0/M0 = 4 7 M c M p /% 



(8) 



This energy is directly related to the strength of the in- 
terparticle interactions. 
Solving Eq. o we obtain 



it = u c + In ( 1 



(9) 



where u c — U c /T is the reduced effective energy bar- 
rier at t = 0, the time when the relaxation starts. The 
normalization time to is given by 



,'U C \ T /Uo-© 
t„^exp[-j=— expl^^ 



(10) 



Now, using Eqs. I and I we get the time evolution of the 
magnetic moment: 



M (t) = M c 1 



T f t 

e" 1 ( 1 + T a 



("J 



Normalized relaxation rate R = \dM/dhx [t)\ /M c is 
given by: 



R = 



T t 

eto + t 



(12) 



As we will see below, experiment shows that to < 1 sec. 
In our measurements typical time window At w 100 sec, 
therefore we can assume t 3> to and Eq. |l2| predicts 
that the relaxation rate saturates at R = T/Q. Thus, 
measurements of the normalized relaxation rate can pro- 
vide direct estimate of the energy scale governing the 
relaxation process. 



C. Irreversibility temperature 

The irreversibility temperature Ti rr of the assembly 
of magnetic nanoparticles is defined by the condition 
M (At, T irr ) = AM. Here AM is the smallest measured 
magnetic moment and At is the time window of the ex- 
periment. Using Eq. ^ we obtain: 



T 



e 







1 - AM/M C 
In (1 + At /t ) ~ In (1 + At /to) 



(13) 



thus we can estimate the characteristic time to from 
measurements of Ti rr , because the energy can be de- 
termined separately from the measurements of the re- 
laxation rate R w T/Q. On the other hand the irre- 
versibility line T irr (H) gives the field dependence of 0. 
The latter may be obtained also from R (H) measure- 
ments. Thus measurements of Ti rr (H) and R (H) in 
different samples provide a verification of our model on 
self-consistency. 

It is interesting to note that the expression for Ti rr , Eq. 
[l3] , is typical for the blocking temperature of individual 
non-interacting particles, which is obtained from Eq. ^[ 



In (At/t*) 



(14) 



where t* = 1/ (Ahi \M C /AM\) is the characteristic time 
and Uo is given by Eq. |]. Energy Uo is proportional to 
KV for non-interacting nanoparticles, but it is reduced 
by a term proportional to 7 due to interparticle interac- 
tions. This is in agreement with previous works where 
"static" modifications of the barrier for relaxation were 
consideredtllOEj. Irreversibility temperature, Eq. |l4| 
approaches when AM — > 0, and so does T® rr . This 
reveals an important difference in the physics of the ir- 
reversibility line in interacting and non-interacting par- 
ticles. In the former, there is a true irreversibility in 



5 



the limit AM — > associated with freezing of magnetic 
moments due inter-particle interactions. In the case of 
non-interacting nanoparticles, the apparent irreversibil- 
ity is due to experimental limitations (finite sensitivity, 
e.g., AM). It is important to stress that this is true only 
on a macroscopic time scale At 3> t* , such as relaxation 
or M (T) measurements. If, however, At < t* is real- 
ized, for example in Mossbauer measurements, one may 
detect. ±he irreversibility temperature according to Eq. 

We also note that in any case Tj rr is a dynamic 
crossover from reversible to irreversible state and is de- 
fined for a particular experimental time window At. 

In the following section we compare our experimental 
observations with the model developed above. 



of the irreversible behavior is the largest for a crystalline 
sample. The amorphous sample containing acicular par- 
ticles occupies the intermediate space and the amorphous 
sample embraces the smallest space in this T ~ H phase 
diagram. Such behavior is naturally explained in terms 
of a strength of inter-particle interactions, which are the 
smallest in the case of a regular amorphous sample, inter- 
mediate for an acicular amorphous sample (due to shape 
anisotropy) and the largest for a crystalline sample due to 
crystalline anisotropy. Also, highest irreversibility tem- 
perature of crystalline sample is understood on the basis 
of its largest particle size. 

V. SUMMARY AND CONCLUSIONS 



IV. DISCUSSION 

Our phenomenological model provides a description of 
the irreversible magnetic behavior in the assembly of fer- 
romagnetic nanoparticles. In particular, the model pre- 
dicts the time-logarithmic deca y o f the magnetization, 
see Eq. |ll|. Also, Eqs. 12 and 13 relate both the irre- 
versibility line and the relaxation rate to a single param- 
eter 6 = A'yMcMp/n. 

The magnetic relaxation data of Fig. 5 reveal, indeed, 
time-logarithmic relaxation. Fitting these data to Eq. 
[n] yields the parameters of M c and 9. In Fig. 14 we 
plot the derived energy as a function of temperature 
for the three samples and find that O is the largest for 
a crystalline sample, intermediate for an acicular amor- 
phous and the lowest for a regular amorphous sample. 
The straightforward explanation is that in a crystalline 
sample both 7 and M p are the largest; in an acicular 
amorphous sample M p is of the same order as in regular 
amorphous, but 7 is much larger due to shape anisotropy. 

Similarly, we derive the magnetic field dependence of 
from the data of Fig. 9 and plot (H) for three samples 
in Fig. 15. We note that M p and 7 should not depend 
on magnetic field. It is therefore expected that the field 
dependence of is determined by the field dependence 
of M c w M s — M rev (H), which decreases with field, Fig. 
9. Figure 15 shows the agreement with this observation. 
The weak increase of with temperature, Fig. 14, may 
be related to some non-linear dependence of barrier U on 
the magnetic moment. 

Independent estimations of are derived from Ti rr 
of Fig. 7 using Eq. [ll| Comparing Fig. 7 and Fig. 
15 we get (for three samples) <d/T irr « 4 — 6. Thus, 
t ~ 0.05 — 0.5 sec. Note that these values of to are 
much-larger than the "microscopic" values predicted by 
NeelED, simply because they reflect collective behavior of 
the whole assembly controlled by the effective barrier 0, 
see Eq.|l0|, and not a single particle barrier KV. 

Let us now compare the irreversibility lines of different 
samples, Fig. 7. In most parts of this diagram the region 



We presented measurements of irreversible magnetiza- 
tion as a function of temperature, time and magnetic field 
in three types of ferromagnetic nanoparticles: regular 
amorphous, acicular amorphous and crystalline nanopar- 
ticles. The results are interpreted using a developed phe- 
nomenological approach based on the assumption that 
the barrier for magnetic moment reorientation depends 
on the total magnetic moment via dipole fields. This ex- 
plains the time-logarithmic magnetic relaxation governed 
by the energy scale related to interparticle interaction. 
Values of found from measurements of the irreversibil- 
ity line and the relaxation rate are in perfect agreement, 
implying validity of our model. 
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VI. APPENDIX: EFFECTIVE BARRIER FOR 
MAGNETIC RELAXATION AND EQUATION 
FOR TIME EVOLUTION OF THE MAGNETIC 
MOMENT 

Here we consider in detail the model outlined in the 
text. We assume that magnetic moment M p of any 
given particle i can be in one of the two possible energy 



minima: W{ « -M p Hi cos (ft) or W\ M p HiCOs(9i). 



These minima are separated by the barrier of height 
~KV + M p Hi sin (ft). In the presence of thermal fluctu- 
ations, a particle moment sitting in the minima W{ can 
spontaneously change its direction to the next minima 
W 2 . The energy barrier for such reorientation is 



U{ 2 = KV + M p Hi (sin (ft) + cos (ft)) 



(15) 



The backward rotation is also possible and requires over- 
coming the barrier: 



U l 2l = KV + M p Hi (sin (ft) - cos (ft)) 



(16) 



From this point on one can conduct a self-consistent 
statistical average over angles ipi (-H,-,t, ft) in order to 
evaluate the resulting magnetic moment M of the system. 
On the other hand we may try to simplify the problem 
assuming that the internal field Hi can be represented 
as a simple sum of the external field H and the collincar 
to it dipole field Hd (i.e. independent of ft). If all easy 
axes are randomly distributed the average barrier for flux 
reorientation is then given by 



where k = 12 or 21 denotes particle's moment flipping 
from the minima W\ to the minima W 2 , or backward, 
respectively. Using Eq. fL5l and Eq. 16 



we find 



Uia = KV + AMpHi/ir 
U 2 i = KV 



(17) 
(18) 



Let us now consider a situation where temperature is 
higher than irreversibility temperature and system is at 
thermal equilibrium. The number of particles jumping 
per unit time from one minima to another is proportional 
to Nk exp(— Uk/T). The condition for equilibrium is 



iVie" 



U12 



Noe 



.£21 



Thus 



N 2 = Ni exp 



U12 - U21 
T 



Ni exp 



AMpHj 
irT 



It is clear that the difference n = N\ — N 2 determines 
the resulting magnetic moment of a system. If the total 
number of particles in the system is N difference, n is 



N t - N 2 = N - 



exp 



( 4M P H, \ 
\ * T J 



1 + exp 



( _iM £ H ± \ 



TVtanh 



The total reversible magnetic moment then is 



M rev ss M v n 



MpNtznhf'^Ei 



irT 



1 ttT 



(19) 



(20) 



This formula is similar to the expression for the Ising su- 
perparamagnet rand simply reflects the two-state nature 
of our modelEZIo. The difference is, however, that the 
physical magnetic field is the total (external + dipole) 
field H t . 

Dipole field Hd at any given site depends upon ori- 
entations of the moments of the surrounding particles. 
If those orientations are totally random (minima W\ 
and W 2 are equally occupied) the dipole field is small, 
whereas if all surrounding particles are situated in one of 
the minima the resulting dipole field is maximal. From 
this simple picture, we conclude that the magnitude of a 
dipole field depends upon the total magnetic moment of 
a sample M rev + M, where M rev is given by Eq. ^ and 
M is the irreversible, time dependent contribution to the 
total magnetic moment resulting from the finite relax- 
ation time needed for a system to equilibrate. Therefore, 
we may write Hi = H — 7 (M rev + M). Here 7 is the co- 
efficient accounting for the contribution of dipole fields. 
Now we can obtain the equation for reversible magneti- 
zation from Eqj2fj|: 



tanh ( 



2M„H 



M r 



Ms 



2M P M S 
' 7 ttT 



1 



(21) 



where Ms = M p N. We note that this formula is valid 
at small enough fields 2M p H/ir < T when particle mo- 
ments are almost locked along the easy axes and small 
enough interactions (i.e. H > 7Af ret ,) . The important 
result is that reversible magnetization decreases as the 
inter-particles interaction increases. Interestingly, Eq. |2l| 
provides a good description of the experimental data. 

Thus, the barriers for moment reorientation in Eq. [l7] 
and Eq. [l8] can be re- written as 

Ui = KV + AM P (H - 7 (M rev + M)) /tt (22) 
U 2 = KV (23) 

We shall now consider direct and backward moment 
rotation processes in a non equilibrium state. As above, 
we denote by TVi and N 2 number of moments in energy 
minima 1 and 2, respectively. The total number of parti- 
cles in the system is N = Ni+N 2 . The magnetic moment 
is proportional to the difference n — N\ — N 2 . During 
small time St this difference changes as 



Sn 



Ni exp 



T 



N 2 exp 



T 



St (24) 
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Using simple algebra and above relationships between 
Ni, N 2 , n, and N, we get 



¥ = " eXP 



U1 + U2 
T 



cosh 



U1-U2 
T 



TVsinh 



(25) 



From this we arrive to a non-linear differential equation 
governing process of magnetic relaxation not too close to 
equilibrium: 



dM 

— --Aexp 



U1 + U2 
T 



(M rev + M) cosh 



U1-U2 
T 

(26) 



where A is a constant measured in sec -1 and having 
meaning of attempt frequency. 

Eq. [2^ can be simplified considering magnetic relax- 
ation not too close to equilibrium and retaining our as- 
sumption that anisotropy contribution to the magnetic 
energy is much larger than that of magnetic field (both 
conditions are better satisfied at low fields). In this case, 
Eq. may be approximated in a reduced form: 



dM 
~dt 



-AM c exp(-U/T) 



(27) 



where M c is the total magnetic moment at the beginning 
of relaxation and U is the effective barrier: 



U = 2KV + &M p (H - -/M r 



M 
(28) 



where U = KV + AM p (H — jM rev ) /it and M 
^{^ + H~lM r , 



We reiterate that Eq. |27] is valid only in the case when 
the magnetic anisotropy is large and magnetic moment 
is far from equilibrium. Close to equilibrium, one ought 
to consider Eq. |2(| 
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VII. FIGURE CAPTIONS 



Fig.l X-ray diffraction patterns for (a) regular amor- 
phous, (b) acicular amorphous and (c) crystalline sam- 
ples. 

Fig. 2 Differential scanning colorimetry spectra for (a) 
regular amorphous, (b) acicular amorphous and (c) crys- 
talline samples. 

Fig. 3 Transmission electron micrographs for (a) regular 
amorphous, (b) acicular amorphous and (c) crystalline 
samples. 

Fig. 4 Typical ZFC - FC curves with superimposed re- 
laxation for the amorphous sample measured at 500 G at 
different temperatures. Inset: full-range ZFC-FC curve. 

Fig. 5 Typical relaxation curves measured in the amor- 
phous sample in 500 G at different temperatures. 

Fig. 6 Normalized logarithmic relaxation rate R for 
three types of samples as a function of temperature. In- 
set: R as a function of a reduced temperature T/Ti rr . 

Fig. 7 Irreversibility lines for three types of samples. 

Fig. 8 Typical magnetization loops at T = 5 K (open 
circles) and at T = 100 (solid line). 

Fig. 9 Magnetic relaxation at different values of mag- 
netic field. Vertical lines are the M (t) curves superim- 
posed on a regular magnetization loop measured at the 
same temperature. 

Fig. 10 Normalized logarithmic relaxation rate R for 
three samples as a function of magnetic field at T = 20 K. 

Fig. 11 Schematic snapshots of magnetic moments dis- 
tribution in powder sample at (a) beginning of the relax- 
ation and (b) at latter time 

Fig. 12 Energy profiles after FC in magnetic high field 
(Hi) and after reduction of the magnetic field, whence 
the relaxation starts (H^)- 

Fig. 13 Energy profiles at the beginning and at the lat- 
ter stage of the relaxation. Dots indicate population of 
magnetic moments of the particular energy minima. 

Fig. 14 Temperature dependence of energy extracted 
from the measurements of normalized relaxation rate for 
the three samples. 

Fig. 15 Magnetic field dependence of the energy <d. 

Fig. 16 Scaling of the M (T) FC-ZFC curves with ir- 
reversibility temperature. 
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Table I. Characteristic parameters of the samples 



Sample 


d [nm] 


Surface area [m 2 / g] 


regular amorphous 


~50 


148 


acicular amorphous 


~5 x 50 


164 


crystalline 


~200 


88 1 
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